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ELASTIC  DISPLACEMENT  OF  PRIMARY  WAVES  FROM 
EXPLOSIVE  SOURCES 

By  Joseph  W.  Berg,  Jk.,  and  George  E.  Papageorge 
ABSTRACT 

Equations  derived  from  Blake's  (1952)  model  of  diverging  waves  from  a  point  Bource  were 
programmed  for  computer  analysis.  Variations  of  the  displacement  resulting  from  changes  in 
the  pressure  function,  propagational  velocity,  range,  and  cavity  size  were  investigated.  Re¬ 
sults  of  the  research  indicate  that:  (1)  a  step  pressure  function  used  with  this  model  gives  dis¬ 
placements  that  closely  approximate  the  displacements  measured  at  0.4  km  from  the  Gnome 
nuclear  explosion;  (2)  near  the  source,  long-period  displacements  are  inherent  with  this  model; 
(3)  the  periods  of  the  maximum  Fourier  transform  amplitude  of  the  radiation  field  is  propor¬ 
tional  to  the  equivalent  cavity  radius ;  (4)  the  peak  displacements  scale  to  the  two-thirds  power 
of  charge  size  for  values  of  a  between  80  (0.5  kt)  and  600  m  (276  kt);  and  (6)  between  0.1  and 
3.0  cps,  the  amplitudes  of  given  frequencies  scale  to  the  first  power  of  charge  size  for  values  of 
a  between  145  (3  kt)  and  305  m  (28  kt).  In  general,  Fourier  amplitudes  at  frequencies  below 
the  natural  frequency  of  the  cavity  scale  to  the  first  power  of  charge  size,  and  Fourier  am¬ 
plitudes  at  frequencies  above  the  natural  frequency  of  the  cavity  scale  to  a  fractional  power  of 
charge  size.  It  is  suggested  that  this  may  be  a  good  model  with  which  to  compare  near-source 
observations  of  seismic  phenomena. 


Introduction 

This  investigation  compares  a  theoretical  description  with  results  derived  from 
from  observed  data  for  compressional  seismic  waves  generated  by  an  explosive 
source.  In  this  work,  theoretical  descriptions  of  ground  displacements  are  con¬ 
sidered  in  the  elastic  region  only,  between  the  ranges  of  0.46  and  15.2  kilometers. 

Mathematical  descriptions  of  spherically  diverging  elastic  waves  instigated  by 
pressure  functions  have  been  given  by  Sharpe  (1942),  Duvall  and  Atchison  (1950), 
and  Blake  (1952).  Blake’s  model  was  chosen  for  this  investigation.  In  this  model, 
long-period  ground  displacements  which  decay  exponentially  with  the  pressure, 
and  oscillatory  transients  of  the  radiation  field  are  important  for  the  ranges  con¬ 
sidered.  Very  little  information  regarding  long-period  displacements  is  available, 
but  some  observations  have  been  made  by  Weart  (1962)  and  Werth  and  Herbst 
(1963). 

In  the  following  work,  variations  of  ground  displacement  have  been  considered 
with  changes  in  pressure  functions,  range,  compressional  wave  propagational 
velocity,  and  cavity  size. 


Analyses 

Source 

For  an  infinite  homogeneous,  isotropic  medium,  the  solution  to  the  wave  equa¬ 
tion  is 

¥  -  (A/r)/(r)  (1) 

where  is  a  potential  function  of  displacement,  A  is  a  constant,  r  is  range,  r  « 
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'  —  (r  —  a)/c,  a  is  the  radius  of  the  cavity,  t  is  time,  and  c  is  compressional  wave 
propagational  velocity. 

The  boundary  condition  requires  that  the  radial  stress  in  the  cavity  wail  equals 
the  pressure  within  the  cavity.  This  is  expressed  by 


P(t)  =  -pc8 


(2) 


where:  P(t)  is  the  pressure  function:  p  is  density;  u  is  displacement;  and  a  is  Poisson’s 
ratio,  taken  to  be  0.25  in  this  work.  Blake  obtains  the  following  formal  solution 
of  the  wave  equation  in  terms  of  displacement  potential,  e>,  for  any  given  pressure 
function  P(y) 


<P  = 


(x*  -  jx/k  -  1/A) 


(3) 


where:  x  =  wa/c;  A-  =  |(1  —  <r)(l  —  2<r)-1;  w  is  angular  frequency;  and  y  is  an 
integration  parameter. 

The  pressure  functions 


P(y)  =  Po(e~a,y  -  e~a,y)  for  7  ^  0 


and 


P(y)  =  0 


for  7  <0  (4) 


were  used  to  solve  Equation  (3).  A  pressure  function  of  this  form  allows  choice 
in  shaping  the  pressure  by  varying  the  values  of  co  and  at  between  0  and  «. 

Equation  (3)  was  solved,  and  it  was  then  differentiated  with  respect  to  r  to 
give  the  following  equation  for  displacement. 


u  = 


P oO 


pr[uos  +  (a0 

__  a0  —  ai 
Wo 


)  +  " 


cos  UoT 


sin  wor 


^  ,  (a to  —  ati)ao  +  wo'  , 

cos  wor  H - sin 


in  “bT|J 


Pan 


pr[wos  +  (a0 

ato  —  aj 


Wo 


') 


(5) 


cos  wor 


•„  1  ,  aj _  (ao  —  «i)aro  +  Wi  . 

sin  wor  1  -+•  —  cos  wor  + - sin  wor 

c  Wo c 


}] 


E( illation  (5)  contains  two  types  of  terms:  those  which  have  a  geometrical 
divergence  of  1  /r* ;  and  those  which  have  a  divergence  of  I  r.  Included  in  the  terms 
whose  amplitudes  diminish  as  1  r*  are  those  which  depend  upon  the  exponential 
decay  of  the  forcing  function.  These  terms  represent  long-period  displacements 
and,  in  the  case  of  a  pressure  stepwave,  part  of  these  terms  represent  permanent 
displacement.  The  permanent  displacement  results  from  the  infinite  duration  of 
the  pressure  on  the  wall  of  the  cavity. 

Equation  (5)  was  programmed  for  an  IBM  1020  Computer  to  investigate  the 
effect  of  pressure  shape  and  other  factors  on  the  ground  displacement  of  the  primary 
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waves.  Figure  1  shows  the  approximations  to  the  ground  displacement  measured 
(Weart,  1962)  at  298  meters  from  the  Gnome  nuclear  explosion.  The  approxima¬ 
tions  shown  on  figure  1  were  made  using  Equation  (5),  on  =  0  and  at  =  »,  and 
assuming  that  the  equivalent  radius,  a  (the  radius  beyond  which  the  medium 
behaves  elastically),  for  a  one-kiloton  yield  closely  coupled  to  the  medium  is  100 
meters  (Carpenter,  Savill,  and  Wright,  1962).  It  was  assumed  that  a  scales  as 
IF1'’.  Thus,  the  equivalent  radius  for  Gnome  (3  kt)  was  calculated  to  be  146 


TIME  (»*ca«d>) 

Flo.  1.  Comparison  of  theoretical  ground  displacement  amplitudes  with  that  observed  from 
the  (inome  nuclear  explosion  as  given  by  Weart  (1962). 

meters.  If  it  was  desired,  a  closer  approximation  to  the  observed  data  could  be 
obtained  by  changing  the  values  of  P0 ,  an  ,  a» ,  a,  c,  and  possibly  other  considera¬ 
tions.  From  figure  1,  the  step  pressure  function  [Eq.  (6),  an  **  0,  a*  —  *]  would  be  a 
close  approximation  to  tlie  observed  data.  This  agrees  with  Latter,  LeLevier, 
Martinelli,  and  McMillan  ( 1969). 

Distance  from  Source 

The  displacement  function  at  ranges  of  0.46  km,  fig.  2A,  and  16.2  km,  fig.  2B, 
is  plotted  for  equivalent  cavity  radii  of  146  m  (3  kt  yield)  and  306  m  (28  kt  yield). 
For  these  calculations  the  pressure  was  assumed  to  be  a  step  function.  Long-period 
displacements  are  evident  at  the  0.46  km  range  whereas  they  are  not  as  evident  at 


10 


FREOUENCY  (cps ) 

Fio.  4A.  Fourier  Transform  Amplitudes  versus  frequency  for  different  equivalent  cavity  radii 

at  a  range  of  0.46  km. 


Fio.  4B.  Fourier  Transform  Amplitudes  versus  frequency  for  different  equivalent  cavity  radii 

at  a  range  of  lo.2  km. 
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the  15.2  km  range.  For  computer  calculations,  all  waveforms  were  truncated 
after  1  second.  It  is  believed  that  no  serious  omissions  of  data  resulted  from  this 
operation.  The  Fourier  Transform  amplitudes  for  frequencies  of  0.1,  0.5,  1.0,  and 
3.0  cps  are  plotted  against  distance  for  these  ranges.  Straight  lines  are  used  to 
represent  the  decrease  of  amplitude  with  distance  for  0.5  and  0.1  cps  between  0.3 
and  1.6  km.  '  'he  slopes  of  the  line  segments  are  —  1.9  cm-sec/km.  If  the  displace- 


Fio.  5.  Period  of  peak  amplitude  (see  peaks  marked  T  on  Figures  4A  and  4B)  versus  equivalent 

cavity  radius. 

meat  waveforms  had  not  been  truncated  at  I  second,  the  amplitudes  of  the  lower 
frequencies  (0.1  and  0.5  cps)  would  diminish  by  the  inverse  square  of  the  range, 
and  at  distances  greater  than  15  km  they  would  diminish  as  the  inverse  range 
(see  Equation  (5)j.  The  amplitudes  of  higher  frequencies  (1.0  and  3.0  cps)  diminish 
as  the  inverse  range. 


Amplitudes,  Periods,  and  Site  of  Source 

The  Fourier  Transform  amplitudes  versus  frequency  are  plotted  for  ranges  of 
0.40  km  (fig,  4A)  and  15.2  km  (fig.  4B).  The  long-period  displacements  are  in 
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evidence  for  frequencies  below  about  1  ops;  they  are  considerably  smaller  at  15.2 
km  than  at  0.46  km  from  the  source.  These  long-period  displacements  have  been 
difficult  to  observe,  but  they  have  been  observed  near  nuclear  explosions  (Weart, 
1962.  and  Werth  and  Herbst,  1963).  However,  additional  evidence  for  their  exist¬ 
ence  is  desirable. 


Fio.  (iA.  Peak  displacement  amplitude  versus  equivalent  cavity  radii  for  range  of  O.-Mi  km, 
velocity  of  4.88  km/sec,  and  on’s  of  0,  1,  3,  and  4,  and  o»  —  », 

For  frequencies  between  1  and  10  cps,  the  amplitude  peaks  for  the  radiation 
field  part  of  the  theoretical  source  are  indicated  in  figures  4A  and  4B  by  the  vertical 
lines  marked  "T".  It  is  easy  to  see  the  increase  in  frequency  with  decreasing  « 
(equivalent  cavity  radius).  If  a1  is  proportional  to  W  (charge  sise),  the  shift  in 
period  of  the  peak  would  be  expected  to  be  proportional  to  the  equivalent  cavity 
radius,  a,  as  Latter,  I^eLevier,  Martinelli,  and  McMillan  (1959)  give  the  scaling 
of  characteristic  frequency,  <* ,  as  proportional  to  Figure  5  shows  a  plot  of 
period  of  peak  amplitude  versus  cavity  radius.  These  variables  are  proportional 
at  the  two  distances  considered ;  the  slopes  of  the  straight  lines  vary  with  distance. 
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The  waveform  of  ground  displacement  as  given  by  equation  (5)  depends  upon 
pressure  pulse  attenuation,  ground  velocity,  range,  and  cavity  size.  For  peak 
ground  displacement  amplitude,  A,  the  exponent  n  in  the  relationship  A  <*  W* 
was  investigated  for  a  practical  range  of  the  above  variables.  Figures  6A,  6B,  6C, 
and  6D  show  the  peak  displacement  amplitudes  plotted  against  the  equivalent 


Fig.  OB-  Peak  displacement  amplitude  versus  equivalent  cavity  radii  for  range  of  0.40  km, 
velocity  of  2.44  km/sec,  and  cn’s  of  0,  1,  3,  and  4,  and  a,  —  • . 

cavity  radii  for  values  of  on  between  0  and  4,  and  at  *  «.  Figures  6A  and  6B 
are  computed  for  velocities  of  4.88  and  2.44  km/sec  at  a  range  of  0.46  km;  figures 
6C  and  6D  at  a  range  of  15.2  km.  Computations  to  compare  the  values  of  displace¬ 
ments  for  different  times  along  the  waveform  and  to  read  the  maximum  displace¬ 
ment  were  performed  on  a  computer.  If  a*  is  assumed  proportional  to  W,  values  of 
n  in  the  relationship  A  —  kW9  are  seen  to  vary  between  0.64  and  0.79*  for  values 


*  O’Brien  (1960)  summarised  values  of  the  exponent  a  given  by  different  investigators. 
The  value*  ranged  between  0.73  and  1  JO.  The  largest  charge  sis*  used  to  obtain  these  values 
was  300  lb*.  A  value  of  0.78  for  a  was  given  by  Carder  and  Cloud  (1900)  for  surface  wave*  with 
periods  of  1  to  2  seconds  for  larger  explosions. 
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of  a  between  80  and  600  m.  This  is  not  a  large  variation  of  n.  This  method  of  relat¬ 
ing  the  peak  displacement  amplitude  to  the  size  of  the  source  is  equivalent  to 
using  the  peak  displacement  from  the  first  half  cycle  of  a  displacement-type  seismo¬ 
gram.  These  results  are  comparable  to  those  of  Werth  and  Herbst  (1663). 

It  must  be  realized  that  frequency  shifts  with  different  charge  sizes.  Thus,  it 


Fig.  AC.  Peak  displacement  amplitude  versus  equivalent  cavity  radii  for  range  of  15.2  km, 
velocity  of  4.88  km/sec,  and  oi’s  of  0,  1,  3,  and  4,  and  at  ■  >. 

would  be  expected  that  the  amplitude  of  a  given  frequency  would  scale  differently 
to  charge  size  than  the  peak  displacement  amplitude.  The  Fourier  transform 
amplitudes  versus  the  equivalent  cavity  radii  are  plotted  at  ranges  of  0.46  (fig. 
7A)  and  15.2  km  (fig.  7B).  At  the  0.46  km  range  (fig.  7A),  the  Fourier  amplitudes  of 
frequencies  between  0.1  and  3  cps  scale  to  the  first  power  of  charge  size  (n  1) 
while  those  at  5  cps  scale  to  the  0.78  power.  The  two  different  amplitude  groupings 
in  figure  7  A  are  due  to  the  relative  size  of  the  long-period  displacements  compared 
to  the  amplitudes  of  the  radiation  field  at  this  range.  At  the  range  of  15.2  km  (fig. 
7B),  the  Fourier  amplitudes  between  0.1  and  1  cps  scale  to  the  first  power  of  the 
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charge  size  while  those  at  3  and  5  cps  are  not  linear  in  this  region.  The  natural 
frequencies  of  the  cavities  (fo)  were  calculated  to  be  5  cps  and  3  cps  for  equivalent 
cavity  radii  of  145  and  305  meters,  respectively,  and  for  a  compressional  wave 
velocity  of  4.88  km/sec.  Thus,  for  the  above  range  of  frequencies,  which  are  below 
the  natural  frequencies  of  the  cavities,  the  Fourier  Amplitudes  scale  to  the  first 


Fin.  6D.  Peak  displacement  amplitude  versus  equivalent  cavity  radii  for  range  of  15.2  km, 
velocity  of  2.44  km/sec,  and  ai’s  of  0,  1.  3,  and  4,  and  at  —  ». 

power  ( n  ^  1 )  of  the  charge  size.  This  is  in  accord  with  Latter,  Martinelli,  and 
Teller  ( 1959).  Moreover,  an  analysis  of  Equation  (4)  of  their  work  shows  that  the 
Fourier  Amplitudes  scale  to  the  first  power  of  the  charge  size  for  very  low  fre¬ 
quencies  («  <  wo)  and  to  the  one-third  power  of  the  charge  size  for  high  frequencies 

(»  >  Mo).* 

A  value  of  n,  in  the  relationship  .4  «  IF",  that  is  different  from  0.5  should  not  be 
of  great  concern.  The  explosion  radiates  a  spectrum  of  energy,  and  not  energy  at 
one  particular  frequency  in  the  form  of  a  sine  wave.  The  amplitudes  of  the  different 


In  some  instances  the  exponent  «  may  lie  as  large  as  4jj. 


EQUIVALENT  CAVITY  RADIUS  IN  METERS 


Flo.  7A.  Fourier  Transform  Amplitudes  versus  equivalent  cavity  radii  for  frequencies  of  0.1. 
0.3,  0.5,  1,  3,  and  5  cps  at  range  of  0.4<i  km. 


Fio.  7B.  Fourier  Transform  Amplitudes  versus  equivalent  cavity  radii  for  frequencies  of  0.1, 
0.3,  0.5,  1,  3,  and  5  cps  at  range  of  15.3  km. 
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frequencies  scale  differently  to  charge  size  as  shown  above.  This  was  pointed  out 
by  Peet  (1960). 


Conclusions 

In  general,  the  theoretical  source  investigated  in  this  research  gives  results  com¬ 
patible  with  observations  made  by  other  investigators.  The  results  indicate  that: 
(1)  a  step  pressure  function  used  with  this  model  gives  displacements  that  closely 
approximate  the  displacements  measured  at  0.4  km  from  the  Gnome  nuclear  ex¬ 
plosion;  (2)  near  the  source,  long-period  displacements  are  inherent  with  this  model; 
(3)  the  period  of  the  maximum  Fourier  transform  amplitude  of  the  radiation  field 
is  proportional  to  the  equivalent  cavity  radius;  (4)  the  peak  displacements  scale 
to  the  two-thirds  power  of  the  charge  size  for  values  of  a  between  80  (0.5  kt)  and 
600  m  (215  kt);  and  (5)  between  0.1  and  3.0  cps,  the  amplitudes  of  given  frequencies 
scale  to  the  first  power  of  charge  size  for  values  of  a  between  145  (3  kt)  and  305  m 
(28  kt).  In  general,  Fourier  amplitudes  at  frequencies  below  the  natural  frequency 
of  the  cavity  scale  to  the  first  power  of  charge  size,  and  Fourier  amplitudes  at 
frequencies  above  the  natural  frequency  of  the  cavity  scale  to  a  fractional  power  of 
charge  size. 

The  long-period  displacements  given  by  this  model  for  near-source  ranges  offer  a 
problem  in  that  they  have  been  difficult  to  observe  to  date.  Such  displacements 
could  be  very  important  when  considering  attenuation  of  energy  near  the  source, 
generation  of  surface  waves,  and  other  seismic  phenomena. 
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